#########################################################################################
# The Effect of Conflict-Related Sexual Violence on Civic Engagement
# 
# Figure 3
# 
#########################################################################################

par(mfrow=c(1,3), mar=c(3, 5, 1, 0), cex.axis=.8, cex.lab=.8, cex.main=.8, las=1, oma=c(1,1,1,1), mgp=c(1.5, .5, 0))

plot(c(mean(civic.congo.ind), 
       mean(civic.congo.dir)), 2:1, pch=c(19, 21), cex=1.2, ylim=c(0.5, 2.5), xlim=c(-.2, .5), main="DRC", axes=F, ylab="", xlab="Difference in Probability")

#grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, lty=2, col="darkgrey")
segments(c(quantile(civic.congo.ind, .05), 
           quantile(civic.congo.dir, .05)), 2:1,
         c(quantile(civic.congo.ind, .95), 
           quantile(civic.congo.dir, .95)), 2:1)
axis(2, 2:1, labels=c("List \n  Experiment", "Direct Item"), col="white")
axis(1, col="white")
points(mean(civic.congo.dir), 1, pch=19, cex=.8, col="white")

text(c(mean(civic.congo.ind), 
       mean(civic.congo.dir)), 2:1, round(c(mean(civic.congo.ind), 
                                            mean(civic.congo.dir)),2), pos=3, cex=.8)
box(col="darkgrey")
par(mar=c(3, 2.5, 1, 2.5))

plot(c(mean(civic.liberia.ind), 
       mean(civic.liberia.dir)), 2:1, pch=c(19, 21), cex=1.2, ylim=c(0.5, 2.5), xlim=c(-.2, .5), main="Liberia", axes=F, ylab="", xlab="Difference in Probability")

#grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, lty=2, col="darkgrey")
segments(c(quantile(civic.liberia.ind, .05), 
           quantile(civic.liberia.dir, .05)), 2:1,
         c(quantile(civic.liberia.ind, .95), 
           quantile(civic.liberia.dir, .95)), 2:1)
axis(1, col="white")
points(mean(civic.liberia.dir), 1, pch=19, cex=.8, col="white")

text(c(mean(civic.liberia.ind), 
       mean(civic.liberia.dir)), 2:1, round(c(mean(civic.liberia.ind),
                                              mean(civic.liberia.dir)),2), pos=3, cex=.8)
box(col="darkgrey")
par(mar=c(3, 0, 1, 5))

plot(c(mean(civic.sri.ind), 
       mean(civic.sri.dir)), 2:1, pch=c(19, 21), cex=1.2, ylim=c(0.5, 2.5), xlim=c(-.2, .5), main="Sri Lanka", axes=F, ylab="", xlab="Difference in Probability")

#grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, lty=2, col="darkgrey")
segments(c(quantile(civic.sri.ind, .05), 
           quantile(civic.sri.dir, .05)), 2:1,
         c(quantile(civic.sri.ind, .95), 
           quantile(civic.sri.dir, .95)), 2:1)
axis(1, col="white")
points(mean(civic.sri.dir), 1, pch=19, cex=.8, col="white")

text(c(mean(civic.sri.ind), 
       mean(civic.sri.dir)), 2:1, round(c(mean(civic.sri.ind), 
                                          mean(civic.sri.dir)),2), pos=3, cex=.8)
box(col="darkgrey")
#########################################################################################

